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ABSTRACT 

Gas rich sub-galactic halos with mass M t < 10 75 M Q , while incapable of forming stars due to lack 
of adequate coolants, contain a large fraction of baryonic mass at cosmological reionization. We show 
that the reionization of the universe at z = 10 — 20 has an interesting physical effect on these halos. 
The external radiation held causes a synchronous inward propagation of an ionization front towards each 
halo, resulting in an inward, convergent shock. The resident gas of mass Mj, ~ 10 4 — 10 7 M Q in low spin 
(initial dimcnsionlcss spin parameter A < 0.01) halos with a velocity dispersion a v < 11 km/s would 
be compressed by a factor of ~ 100 in radius and form self-gravitating baryonic systems. Under the 
assumption that such compressed gaseous systems fragment to form stars, the final stellar systems will 
have a size ~ 2 — 40 pc, velocity dispersion ~ 1 — lOkm/s and a total stellar mass of M* ~ 10 3 — 10 6 M Q . 

The characteristics of these proposed systems seem to match the observed properties of halo globular 
clusters. The expected number density is consistent with the observed number density of halo globular 
clusters. The observed mass function of slope ~ —2 at the high mass end is predicted by the model. 
Strong correlation between velocity dispersion and luminosity (or surface brightness) and lack of corre- 
lation between velocity dispersion and size, in agreement with observations, are expected. Metallicity is, 
on average, expected to be low and should not correlate with any other quantities of globular clusters, 
except that a larger dispersion of metallicity among globular clusters is expected for larger galaxies. The 
observed trend of specific frequency with galaxy type may be produced in the model. We suggest that 
these stellar systems are seen as halo globular clusters today. 

Subject headings: cosmology: theory - dark matter - galaxies: formation - galaxies: kinematics and 
dynamics - globular clusters: general 



1. INTRODUCTION 

Zinn (1985) pointed out insightfully that there appears 
to exist two separate populations of globular clusters in 
the Galaxy: one which is spherically distributed, has low 
metallicity with no radial metallicity gradient and low ro- 
tation, the other which is more concentrated towards the 
Galactic center and has high metallicity, a gradient in 
metallicity and high rotation. Such a dichotomy seems 
to be shared by other galaxies (e.g., Ashman & Bird 
1993; Forbes, Grillmair, & Smith 1997; Forbes, Brodie, 
& Huchra 1997; Kundu 1999; Gebhardt & Kissler-Patig 
1999; Barmby et al. 2000; Beasley et al. 2000). The first 
population, the halo globular clusters in both the Galaxy 
and external galaxies, display an additional array of inter- 
esting properties (Harris 1991). First, they are old with 
age > lOGyr. Second, their characteristics are uniform 
across all galaxy types and sizes. Third, their mass func- 
tion resembles a power-law (oc M~ 1J to ~ 2 - ) at the high 
mass end. 

It has been a long standing challenge to explain such 
features. Several important theoretical models includ- 
ing the pre-galactic Jeans mass based model (Peebles & 
Dicke 1968) and the secondary thermal instability based 
model (Fall & Rees 1985) have been put forth but may not 
be without some intriguing difficulties. For example, the 
Jeans mass based model is not naturally expected to pro- 
duce power-law mass function for globular clusters, while 
the thermal instability based model appears to require del- 
icate heating sources (Harris & Pudritz 1994; Meylan & 



Heggie 1997). A fully consistent model for the origin of the 
halo globular clusters remains unavailable at this time. In 
this paper, we suggest that the reionization of the universe 
could trigger a simultaneous formation of sub-galactic sys- 
tems at that epoch, whose characteristics match remark- 
ably well those of the observed halo globular clusters. 

After describing the globular cluster formation theory 
in §2, we discuss some of the expected properties of glob- 
ular clusters which form in the current theory in §3, fol- 
lowed by conclusions in §4. The following cosmological 
model is used throughout: flcDM = 0.26, = 0.04, 
fl t = n CD M + M b = 0.3, A = 0.0, H Q = 65km/sec/Mpc 
and erg = 1.0; the global gas to total mass ratio R = 
n b /ft t = 0.13. 

2. FORMATION OF GLOBULAR CLUSTERS AT 
COSMOLOGICAL REIONIZATION 

In the standard picture of structure formation, the uni- 
verse is thought to be reionized mostly by photons from 
quasars or stellar systems more massive than 10 9 M 
(Haiman, Rees & Loeb 1997; Gnedin & Ostriker 1997). 
Less massive systems, after having produced a trace 
amount of stars (Pop III) at an earlier epoch, can no longer 
form stars due to lack of cooling processes (Haiman, Thoul, 
& Loeb 1996; Haiman, Rees & Loeb 1997; Tegmark et al. 
1997). During the cosmological reionization phase small 
halos with neutral gas sitting idly within (i.e., without in- 
ternal radiating sources) suddenly see a sea of radiation 
approaching. We will show that this radiation field has a 
dramatic dynamic effect on them. While an exact treat- 
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mcnt of this process is quite intricate and would require 
detailed three-dimensional radiation-hydrodynamical sim- 
ulations, we present a simple analysis that we believe cap- 
tures the essence of the process. To make the calculations 
relatively tractable, we assume, without loss of primary 
characteristic features, that the halos are initially spheri- 
cal with the following density profile: 

P( r ) = Pv{^-y 2 for r <r v 

= Pv (^)- 3 for r>r v , (1) 

where p v = W0p c (z) is the density at the virial radius r v 
and p c (z) = p c (0)Ot(l + z ) 3 is the critical density at red- 
shift z [p c (0) is the critical density today]. The virial mass 
of the halo (i.e., the total mass interior to r v ) is 

M v = 4irp v r 3 v (2) 

and the one-dimensional velocity dispersion within the 
virialized region (r < r v ) is 

a v = ^2irGp v rl (3) 

where G is the gravitational constant. We assume that the 
ratio of gas to total mass in the halos under consideration 
is equal to the global baryon to total mass ratio R. The 
total mass within the radial range r v to r > r v (outside 
the virial radius) is 

M(r v — > r) = M v In — . (4) 

The clumpy universe is reionized outside in: low den- 
sity regions are ionized first and higher density regions 
(except those in the most massive, ionizing sources) be- 
come ionized at progressively later times (Miralda-Escude, 
Haehnelt, & Rees 2000; Gnedin 2000b). The reioniza- 
tion phase progresses on a time scale of a Hubble time 
at the redshift in question; the mean (volume weighted) 
radiation field builds up slowly up to a value of approxi- 
mately 10~ 24 erg/cm 2 /hz/sec/sr at Lyman limit. An ac- 
celerated phase follows, when about 90% of the baryons 
have been ionized and a sudden jump in the amplitude 
of the mean radiation field intensity at Lyman limit to 
10~ 22 — 10~ 21 erg/cm 2 /hz/sec/sr occurs within a redshift 
interval of a fraction of unity and the reionization pro- 
cess is said to be complete (Miralda-Escude et al. 2000; 
Gnedin 2000b). The mean free path of ionizing photons 
approaches the horizon size at a much later time. As Fig- 
ure 3 (below) will show that the gas mass contained in the 
halos of interest here (including gas surrounding the halos 
that will be compressed) is of order 50% of total gas, so 
we are concerned with the phase when a large fraction of 
the gas is not yet ionized and the radiation field is still 
increasing very slowly (see Figure 2a of Gnedin 2000b). 

Let us now examine the global ionization front (I- front), 
which propagates from low density regions towards halos. 
The speed of the external global I-front towards a halo is 

JhR \r v J 

at r > r v , where J v = j°° &irj v dv/hv is the num- 
ber of ionizing photons per cm 2 sec and j„ is the 



intensity of the meta-galactic radiation field; where 
fn = 0.76 is the hydrogen mass fraction of baryons; 
trip is the proton mass. Adopting a powerlaw j„ = 
Jll. 21 10 21 (i'/Vo) - ^ erg/cm 2 /hz/sec/sr one obtains J v = 
1.9 x 10 6 /3- 1 j LL<2 i cm" 2 sec" 1 . Using (3 = 2 and R = 0.13 
we obtain 

where z ri is the redshift when this occurs. Note that for a 
softer stellar spectrum of j„ , (3 ~ 4 and Vj would be lower. 

The I-front is accompanied by a shock front. In or- 
der for the I-front to drive a strong compressive shock 
inward, the shock velocity has to exceed the velocity of 
the ionization front (D-type), which decreases with de- 
creasing halo-centric radius. Otherwise, the neutral gas 
gets ionized and raised to a high pressure before the gas 
is reached by the shock front (see Spitzer 1978 for an in- 
troduction on the subject and Bertoldi & McKee 1990 for 
a detailed treatment in the context of interstellar clouds). 
We assume that this condition sets the initial radius of 
the inward shock, r,: at the ionization front and the 
shock have the same velocity. In the present case the 
shock velocity will instantaneously settle to a speed of 
v shk = A/3^kT„/m p = 12(T H /10 4 fQ 1 / 2 km/s at the 
start of the shock due to the ram pressure (see equation 7 
below), where T r i is the temperature of photo- ionized gas 
and k is the Boltzmann constant. We will return to this 
when we discuss gas cooling in the shell later in this sec- 
tion. Since v s hk > 15 km/s, this initial radius is outside 
the virial radius (see equation 6) (note that, for a more 
realistic density profile where the density slope near the 
viral radius is about —2.4 instead of —3 assumed here, the 
initial shock radius would be still larger). 

We make some simplifying assumptions to treat such a 
shock. We assume that the shock is strong, which is justi- 
fied since neutral gas has a very low temperature, and the 
material that is swept up along its way is accumulated in 
a thin shell obeying the momentum equation 

d{M s v s ) GM{<r)M s 2 

-77 = 5 + ^r z (p ext - p mt ), (7) 

at r z 

where M s is the gas mass within the shell; v s is the veloc- 
ity of the shell; M(< r) is the total mass interior to radius 
r; the first term and the two terms within the parentheses 
on the right hand side are due to gravity, external and in- 
ternal pressure, respectively. In addition, we assume that 
dark matter does not respond to the motion of the gas, 
which is unlikely to be in serious error for our purpose. 
In any case, our assumption is conservative in the sense 
that inclusion of dark matter response to the gas collapse 
would enhance the collapse of halo gas. 

It is instructive to qualitatively consider the calculation 
in two separate regions: r v < r < Ti and < r < r v . At 
r v < r < ri we assume that the gas is initially cold (i.e., 
p in t may be ignored) and p ext = Rp v (ri/r v y 3 kT ri /m p 
is constant. The shell will be subject to external pres- 
sure and gravity, counter-balanced by ram pressure, and 
will obtain a certain velocity when reaching r — r v . At 
< r < r v , if the shell travels supersonically, the down- 
stream pressure is undisturbed. Since initially the gas in 
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the virialized region at < r < r v is in hydrostatic equi- 
librium, the region interior to the shell at any moment 
will remain in hydrostatic equilibrium. In other words, 
for a shell traveling supersonically, the gravitational force 
and interior pressure force just cancel out. In the absence 
of the external pressure and the self-gravity of the shell 
(which is small until the shell becomes self-gravitating), 
the shell will cruise inward (as long as it remains super- 
sonic), except that it slows down due to the swept-up gas. 

Quantitatively, we integrate the momentum equation 
(7) numerically from to r — 0, with the initial condi- 
tion v s (ri) — and M s = 0. We define r/ in the following 
equation and require 77 to be greater than 1 (momentum 
condition) in order for the shell to reach the center of the 
halo: 

mm[v s (r) 



V 



> 1, 



(8) 



where min[w s (r)] is the minimum velocity of the shell in 
the radial range [0, r v ], and a v , the velocity dispersion, is 
used as the isothermal sound speed of the gas at r < r v . 
Note that if equation (8) is not met and the shell slows 
down to a subsonic speed at some radius (likely at a large 
radius since most of the mass is at large radii), the inte- 
rior gas will adjust its pressure causing the interior gas to 
be adiabatically compressed, which will quickly build up 
the internal pressure and cause the shell to slow down and 
eventually to stop. 
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Fig. 1. — shows r) (equation 8; solid curve) as a function of halo 
velocity dispersion a v , with a temperature of the photoionized gas 



of T ri = 2 X 10 4 Kelvin and j LL>2 i = 2 X 10^ at z r 



15. The 



long dashed curve is the total collapsed baryonic mass M gaa (see the 
right vertical axis) and the dot-dashed curve is Mqc = 0.2M gas , 
where 0.2 is an adopted star formation efficiency. 

Figure 1 shows r\ (equation 8) as a function of a v for 
T ri = 2 x 10 4 Kelvin and j LL>2 i = 2 x 1CT 3 at z n = 15. The 
adopted value of photoionized gas temperature is conser- 
vative [Miralda-Escude & Rees (1994) note that the post 
ionization temperature may be as high as 5 x 10 4 Kelvin]. 
We see that halos with velocity dispersion u v < 11 km/sec 
could collapse momentum-wise. The dashed curve is the 
total collapsed baryonic mass M gas (see the right vertical 
axis) and the dot-dashed curve is Mac — 0.2M ffas , where 



0.2 is a putative value for star formation efficiency (see be- 
low for further discussion). We see that compact, baryonic 
systems with mass M gas < 10 6 M© may form. 
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Fig. 2. — shows the initial virial radius (solid curve), the radius of 
the compressed gas cloud (short dashed curve) and the final radius 
of stellar system (globular cluster) (long dashed curve), all in proper 
units, as a function of velocity dispersion at z r i = 15. 

Let us ask how much compression is required to produce 
globular cluster like systems. Figure 2 shows the initial 
virial radius (solid curve) and the radius of the compressed 
gas cloud whose velocity dispersion is a cm — <r v /0.2 (short 
dashed curve) as a function of halo velocity dispersion at 
z r i — 15. We see that the gas clouds need to be compressed 
by a factor of ~ 100. 

A compressed gas cloud will be flattened (due to angular 
momentum) and shocked, cool, presumably fragment and 
form stars. However, we know little about high redshift 
star formation in such systems observationally or theoreti- 
cally, and simply assume that the star formation efficiency 
is 20%. From an analysis of Local Group dwarf galax- 
ies Gnedin (2000a) concludes that star formation at high 
redshift is about 4%, if it was continuous, but could be 
significantly higher, if most of the star formation occurred 
well before lOGyr ago. Therefore, our choice of star forma- 
tion efficiency seems plausible but uncertainty may be as 
large as a factor of two to five. Note that the mean bary- 
onic number density in the compressed cloud is of order 
1000 cm -3 , with the actual final density in a flattened disk 
probably even higher, matching that of molecular cloud 
cores (Friberg & Hjalmarson 1990; Goldsmith 1987). With 
such a large compression the dark matter would become a 
small fraction (of order 1% for a more realistic halo density 
profile) of the final total mass. 

The stars formed will likely blow away most of the re- 
maining gas by supernova explosions (Dekel & Silk 1986; 
Mac Low & Ferrara 1999). After the remaining gas is lost 
as a wind, the cloud will expand by approximately a factor 
of ~ 1/0.2 and result in a final stellar system with a veloc- 
ity dispersion of uqc = 0.2 x a cm = a v and a size shown 
as the long dashed curve in Figure 2. We see that the 
typical size of such a stellar system is between 2 — 40 pc. 
Both the assumption that the final stellar system has the 
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same velocity dispersion as the initial host halo and the 
fudge factor 1/0.2 in setting a cm are rather uncertain, al- 
though these two parameters are degenerate and can be 
condensed into a single parameter. But small variations 
would not change the results qualitatively. Nevertheless, 
detailed simulations will be necessary to definitively quan- 
tify them. 

It is yet unclear whether the required gas compression 
indicated in Figure 2 is achievable in real halos, even if the 
momentum condition rj > 1 (equation 8) is met. The crit- 
ical issue here is that a finite angular momentum of the 
halo gas before compression (Peebles 1969; White 1984) 
could set a maximum compression factor. Let us denote 
the initial spin parameter of the halo (both gas and dark 

matter) as A; = J i \E i \ 1 / 2 G- 1 Mr^ 2 (where J h E { and 
Mi are the initial angular momentum, total energy and 
total mass, respectively), then the final spin parameter of 
the gas cloud is (assuming no loss of angular momentum) 
A/ « \iR^ x a cm j a v Rj 40 Aj (for the short dashed curve 
in Figure 2). Therefore, a typical halo with Aj ~ 0.05 
(Barnes & Efstathiou 1987; Ueda et al. 1994; Steinmetz & 
Bartelmann 1995; Cole & Lacey 1996) would give A/ <~ 2. 
Clearly, only halos with spin parameter values significant 
lower than the typical value can be adequately compressed 
to form the globular cluster like systems. Are there enough 
number of such low spin halos at z ~ 15? 
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Fig. 3. — shows the mass function of halos at z = 15 (solid 
curve) using the Press-Schcchtcr (1976) formula with the follow- 
ing cosmological model parameters: tlcDM = 0.26, Qb = 0.04, 
A = 0.0, Ho = 65km/sec/Mpc and <rg = 1.0. The vertical line 
indicates the halo mass with rj = 1 as shown in Figure 1. The dot- 
ted curve is the cumulative mass fraction in halos, /(> M), and 
the dashed curve indicates the slope of a power-law mass function, 
n(M)dM oc M- 2 dM. 

Figure 3 shows the cumulative fraction of total mass in 
halos at z = 15 (solid curve) using the Press-Schechter 
(1976) formula. The region left of the vertical (long 
dashed) line has rj > 1 (equation 8). We see that the 
(total) mass fraction of halos that meet the condition in- 
dicated by equation (8) is approximately 15%. As numer- 
ical simulations have shown that the distribution of the 
Aj is quite broad with the 80% range being ~ 0.01 — 0.10 



(Ueda et al. 1994; Steinmetz & Bartelmann 1995; Cole & 
Lacey 1996); there are 10% of halos having Ai < 0.01, 
which would give A/ < 0.4 in this model [when being 
compressed by a factor of ~ 100 to r cm (short dashed 
curve) in Figure 2]. If 10% of the halos form globular 
clusters, the resultant density parameter is £lGC,comp ~ 
0.2 x 1.7 x x 15% x 10% w 2 x 10~ 4 (where the 
five terms from left to right are the star formation effi- 
ciency, the ratio of the total swept-up gas mass over the 
gas mass within the virial radius, the baryon density, the 
mass fraction of halos in question and the fraction of low 
spin halos that collapse). The observed ratio of globu- 
lar cluster mass over total Galactic baryonic mass today 
is approximately 0.25% (Harris & Racine 1979), giving 
fico&s = 0.0025fi* = 1.5 x 10~ 5 , where fi* = 0.006 is 
the stellar mass density (e.g., Gnedin & Ostriker 1992). 
Evidently, there are enough low spin halos in the indi- 
cated mass range that could account for the observed halo 
globular clusters. Since the slope of the predicted mass 
function is also in agreement with observations, one can 
infer that the expected number density of globular clusters 
is consistent with observations. The issue of exactly what 
the threshold of the spin parameter is in order to form 
a globular cluster system can only be settled by detailed 
simulations. Complications such as shell instability may 
play a role in this regard. It may be that the initial floe is 
significantly larger than the presently observed value, be- 
cause a large fraction (perhaps relatively faster spinning) 
of the initial globular clusters have largely evaporated or 
disrupted in time and now are a part of the halo stellar 
population. 

The assumption of a thin shell needs further examina- 
tion. With the usual shock jump conditions and equa- 
tion (7), we find that at the start of the shock propaga- 
tion the shell velocity settles to v s = kT r i /m p and the 
shock front travels ahead of the shell with a velocity of 
of v s hk = 4u s /3. The postshock gas temperature would 
be T s = l/3T r j = 6667Kelvin and the postshock pressure 
would be 4p ext /3. Therefore, while the shell is bounded 
at the inner side by the ram pressure, the outer boundary 
will slightly expand initially. In order to keep the shell 
thin the shell gas needs to be cooled. We compute the 
gas cooling including all relevant processes (Haiman et al. 
1996) with the assumption that the gas has been subject 
to the radiation for a period of tnubb (the Hubble time 
at z r i) and the density is Ap v . For the gas in question 
the primary cooling process is molecular hydrogen cool- 
ing and the primary heating process is photo-heating of 
hydrogen atoms. Figure 4 shows the ratio of the cool- 
ing time to the time scale of the system in question t sys 
as a function of optical depth at the Lyman limit, where 
t sys = lkpc/lOkm/scc Rj 1 x 10 8 yr. The incident radia- 
tion field has Jll,21 = 2x 10~ 3 with an index —2 and the 
postshock gas temperature is assumed to be 5 x 10 3 Kelvin. 
We see that gas with a moderate optical depth (shielding) 
of r ~ 4 — 10 can cool efficiently (gas at low end of the 
optical depth is heated up by photo-heating, while gas at 
the high end of the optical depth can not cool efficiently 
due to a low abundance of molecular hydrogen). For a 
shell of neutral gas we can relate its mass to the optical 
depth as M(r) = 3.4x 10 4 r(r/lkpc) 2 M Q . Apparently, ad- 
equate shielding may be achieved for the gas and gas may 
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be able to cool efficiently. Note that the gas cooling time 
remains approximately constant, assuming it cools isobar- 
ically, until it reaches a few hundred Kelvin. In addition, 
self-gravity of the shell will provide some confinement to 
keep the shell thin. More definitive answers can not be 
given without a detailed radiation-hydrodynamic simula- 
tion. 
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Fig. 4. — shows the ratio of the cooling time to the time scale of 
the system in question as a function of optical depth at the Lyman 
limit. The incident radiation field has Jll,21 = 2 x 10 — 3 with an 
index —2. A negative ratio means net heating. The hatched region 
indicates the range in optical depth where cooling is efficient. 

Finally, we give a dimensional analysis for the mass 
threshold of compressed gas clouds. We use the follow- 
ing line of reasoning. The upper mass threshold of the 
cloud should be related to the temperature of the exter- 
nal ionized gas by kT r i = GM v m p /r v , where the virial 
radius r v is related to M v by 5^-r 3 (\ + z r i)Vt t p c (0) — M v 
(5 ~ 200 is the virial overdensity) . The temperature of 
the external ionized gas should be related to hydrogen 
ionization potential kT ri — ^m e c 2 a 2 , where m e is elec- 
tron mass, c is the speed of light and a is the fine struc- 
ture constant. Combining these relations and the Fried- 
mman equation we obtain M gas = 25 f 27r (^r) 3//2 = 

3 x 10 8 i?(i±f^)- 3 / 2 ^7 1/2 /i- 1 M Q , where H is the Hub- 
ble constant at z„ and R is the global gas to total mass 
ratio. It is seen that the upper mass threshold is around 
~ 10 7 - 5 M Q , consistent with the detailed calculations given 
above. 

3. CHARACTERISTICS OF THE GLOBULAR CLUSTERS 

The dotted curve in Figure 3 is the computed mass func- 
tion of halos, n(M) (in units of h 3 Mpc~ 3 ), using the Press- 
Schcchtcr formula and the dashed curve a power-law mass 
function, n(M)dM oc M~ 2 dM. A globular cluster mass 
function of n(M)dM = M~ 2 dM seems to be borne out 
naturally in this picture, in excellent agreement with ob- 
servations for globular clusters more massive than 10 5 Mq 
(Surdin 1979; Racine 1980; Richtler 1992). The observed 
shallower slope at lower mass may be due to disruptive 
effects during subsequent dynamical evolution of globu- 



lar clusters, which tend to work against low mass sys- 
tems (Chernoff, Kochanek, & Shapiro 1986; Chernoff & 
Shapiro 1987; Spitzer 1987; Aguilar, Hut, & Ostriker 1988; 
Chernoff & Weinberg 1990; Harris 1991; Goodman 1993; 
Gncdin & Ostriker 1997). We note, however, there may 
be a cutoff in the velocity dispersion of globular clusters 
at oqc ~ 2km/s, corresponding to a cm <~ lOkm/s, below 
which cooling and thus star formation are less efficient at 
the absence of atomic cooling. 

It is likely that a significant fraction of the initial glob- 
ular clusters may possess a significant amount of angular 
momentum and their initial shapes may be significantly 
flattened. While much more work is needed to quantify 
the dynamical evolution of rotating many-body stellar sys- 
tems (e.g., Einsel & Spurzem 1999), Akiyama & Sugimoto 
(1989) show that initially rotating stellar systems tend to 
become rounder with time. Although it is not possible 
at present to make useful comparisons with observations, 
some initial rotation of the globular clusters might be pre- 
ferred by extant observations (e.g., White & Shawl 1987), 

Simulations of Ueda et al. (1994) (also Cole & Lacey 
1996) show a correlation between average spin parame- 
ter \ and halo mass, Ai oc M~°- 2±01 , for massive (galaxy 
size) systems that they simulate. If we assume that for_the 
halos in question a similar correlation exists, such as \ oc 

M v 1 , and that gas clouds of all halos are compressed 
to a state of comparable final spin parameter, then the 
average final cloud radius fee oc r v X 2 oc r v My X ^ z oc M° 

1 /3 

(because r v oc M v ), i.e., the average final radius of the 
globular clusters is independent of the initial mass of the 
halo, in agreement with the observed lack of correlation be- 
tween the half-light radius and the luminosity (Djorgovski 
& Meylan 1994). Of course, since there is a range in the 
initial spin parameter values, there will be a correspond- 
ing spread in the sizes of globular clusters, as observed. 
Since the final velocity dispersion of a globular cluster is 
cgc fx {My/fh) 1 / 2 , the independence of fh of M v dictates 
ace fx My^ 2 oc L 1 / 2 , in agreement with the observed cor- 
relation ogc « L°- 6±015 (as well as the observed correla- 
tion with surface brightness ace fx /°- 45±0 05 ) (Djorgovski 
& Meylan 1994). A prediction that may be made is that 
the logarithmic dispersion in ace in the ace — L cor- 
relation should be about half the size of the logarithmic 
dispersion in the size of globular clusters, which seems to 
be in agreement with observations (Djorgovski & Meylan 
1994). 

There appears to exist a tendency for halos with rela- 
tively large A to reside preferentially in low density regions 
in simulations (e.g., Ueda et al. 1994). This effect would 
conceivably produce the observed correlation of the spe- 
cific frequency of globular clusters with galaxy type in that 
early type galaxies have a high number of globular clusters 
per unit galactic luminosity than late type galaxies (Har- 
ris 1991). This arises because early type galaxies originate 
from initial higher density peaks and would have (on aver- 
age) lower spins and thus higher specific frequencies. Quite 
intriguing is that this trend might have some bearing on 
the formation of spiral and elliptical galaxies and Hubble 
sequence. But within a similar galaxy type, the number of 
globular clusters in a galaxy may simply be proportional 
to the mass or luminosity of the galaxy, as observed (Hanes 
1977; Harris & Racine 1979), since the mass of a galaxy 
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basically indicates the size of the initial (comoving) region 
from which it has collected matter hence proportionally 
the number of globular clusters. More certain is that the 
range of properties of globular clusters such as size, lu- 
minosity and velocity dispersion should be quite uniform 
across all galaxies, independent of size, type, age, etc, as 
observed (Harris 1991). This is due to the fact that the 
formation of globular clusters in the present model has lit- 
tle correlation with later large-scale density fluctuations 
that form galaxies (except the correlation between spin 
and overdensity mentioned above). 

Chemical enrichment by Pop III stars, which is com- 
puted to be ~ 1(T 3 - 5 Z Q (Ostriker & Gnedin 1996), is 
much lower than the observed stellar metallicity of metal 
poor halo globular clusters of < lCT 1 - 2 Z© (Zinn 1985). It 
thus appears that some low level of self-enrichment has to 
occur in order to account for the observed metallicity of 
globular clusters. While a detailed treatment of star for- 
mation in the proposed globular clusters is not possible, 
the two-generation star formation scenario (Cayrel 1986; 
Brown et al. 1991,1995; Zhang & Ma 1993; Parmentier et 
al. 1999) could work in this picture. We note that the es- 
cape velocity of a compressed cloud (where star formation 
occurs) is larger than its corresponding globular cluster by 
a factor of 5 (this is rather uncertain, though). Therefore, 
confinement of the first generation of supernovae (Dopita 
& Smith 1986) may be more readily accommodated than 
in other scenarios. Then, the second generation of stars 
formed in a massive starburst (with the assumed formation 
efficiency of 20%) completely blows away the remaining 
gas, except for perhaps the most massive globular clusters 
(Dckel & Silk 1986; Mac Low & Ferrara 1999), for which 
further sclf-cnrichmcnt may take place. This is in accord 
with the observed internal chemical homogeneity of most 
globular clusters (cf. Larson 1988), with the exception 
of the most massive systems such as ui Ccn where some 
chemical inhomogencities are observed (Cohen 1981). 

While halo globular clusters depicted here should be, on 
average, metal poor, as observed (Zinn 1985), complex star 
formation histories quite likely introduce a large dispersion 
in metallicity among globular clusters. Thus it may be ex- 
pected that larger galaxies that formed by collecting mat- 
ter in a larger (comoving) region are expected to display 
a larger variance in metallicity, in agreement with obser- 
vations (Harris 1991). Since the metallicity distribution 
is asymmetric and bounded by Z = at the low end, a 
large variance gives a larger mean, interestingly consistent 
with observations (Harris 1991). On the other hand, since 
metallicity does not play a noticeable role in the forma- 
tion of individual globular clusters (i.e., dynamic collapse 
of the gas clouds) in this model, it should be expected to 
show little correlation with any other quantities of globu- 
lar clusters, such as luminosity, size, velocity dispersion, in 
agreement with observations (Djorgovski & Meylan 1994). 

If the universe went through a second, separate reion- 
ization phase of singly ionized helium at lower redshift 
(z ~ 3), as recent observations hint (e.g., Reimers et al. 
1997), a similar radiation I-front will sweep through ha- 
los of neutral gas. However, due to much lower density 
and lower abundance of helium (compared to hydrogen) at 
lower redshift, the speed of the I-front is likely to be larger 
than the shock speed until well within the virial radius; 
i.e, the I-front will remain as an R-type and no significant 



compression will occur. Therefore, our model indicates 
that the age spread of metal poor globular clusters should 
be less than 1 Gyr, consistent with recent observations 
(Rosenberg et al. 1999). 

Finally, it seems inevitable that a fraction of the orig- 
inal globular clusters may be in the intergalactic space 
where no galaxies were formed to collect them. These in- 
tergalactic globular clusters are possibly at a somewhat 
younger dynamic state than Galactic halo globular clus- 
ters due to lack of external dynamical effects. The exact 
number density of intergalactic globular clusters depends 
sensitively on the angular momentum distribution of halos 
in low density regions and requires more detailed calcula- 
tions for quantification. No doubt it will be important 
to search for them with the next generation of large tele- 
scopes. 

4. CONCLUSIONS 

Radiation from stars in galaxies with mass M t > 10 9 M© 
or quasars are thought to be mostly likely responsible for 
reionizing the universe at redshift z ~ 6 — 20. However, a 
large fraction of mass resides in density concentrations on 
smaller mass scales, i.e., gas rich sub-galactic halos with 
mass M t < 10 8 M Q . These smaller systems can not cool 
to form stars without an external trigger. Reionization 
of the universe has an interesting dynamic effect on these 
halos. The rise of the external radiation field causes a 
synchronous propagation of an inward, convergent shock 
towards each halo. We show that, for halos with a ve- 
locity dispersion a v < 11 km/s, the resident gas of mass 
Mb ~ 10 4 — 10 7 M Q will be compressed. For a signifi- 
cant fraction of these halos that have relatively low ini- 
tial angular momentum (A < 0.01), a compression of a 
factor of ~ 100 in radius is possible and compact self- 
gravitating baryonic systems would form. Such gaseous 
systems fragment to form stars and become stellar sys- 
tems of size ^2 — 40 pc, velocity dispersion <~ 1 — lOkm/s 
and a total stellar mass of M s ~ 10 3 — 10 6 M©. 

Their expected properties all seem to be in agreement 
with the observed salient features of halo globular clusters, 
as discussed in the previous section. Angular momentum 
distribution of small halos at high redshift plays a critical 
role in determining many of the properties of globular clus- 
ters. The expected number density of halo globular clus- 
ters is consistent with what is observed. The mass func- 
tion of slope ~ — 2 at the high mass end is predicted by the 
model. Strong correlation between velocity dispersion and 
luminosity (or surface brightness) and lack of correlation 
between velocity dispersion and size are expected. Globu- 
lar clusters should display similar properties, regardless of 
the host galaxy type, age, size or luminosity. Metallicity 
is, on average, expected to be low and should not correlate 
with any other quantities of globular clusters, except that 
a larger dispersion of metallicity among globular clusters is 
expected for larger galaxies. A trend of specific frequency 
with galaxy type may be produced in the model. We pro- 
pose that these stellar systems represent the initial state 
of the presently observed halo globular clusters. 

It will be important to systematically search for high 
redshift as well as intergalactic globular clusters, because 
they will provide tests of the proposed model and could po- 
tentially shed light on the theory of universe reionization 
and the general formation picture of galaxies. Also im- 



7 



portant is to make detailed radiation-hydrodynamic sim- 
ulations with high resolution and superior shock resolving 
power to study the reionization process and follow the dy- 
namics of gas collapse. Both are formidable challenges. 
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